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ABSTRACT 
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Inelastic electron scattering from 
excitation energy revealed giant multipole resonances. Angular distri- 
bution studies with incident 90 MeV electrons show evidence of previ- 
ously unreported giant monopole states. Reduced matrix elements have 
been extracted and multipolarity assignments have been made. Seven 
states have been observed at excitations of 7.3 (undetermined magnetic), 


9.2 (EO), 10.8 (£2), 14.0 (£1), 18.0 (undetermined), 22.5 (E2), and 
33.5 MeV (EO or EZ). 
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I. INTRODUCTION 


Giant resonance investigations of heavy nuclei have been actively 
pursued at the Naval Postgraduate School electron linear accelerator 
since 1973. Initial studies by Warshawsky and Webber [Ref.1] demonstrated 
the feasabiltiy of measuring inelastic electron spectra in the continuum 
region of 19744, They adapted analysis techniques developed at Darmstadt 
and measured cross sections for several known resonances. Warshawsky 
and Webber further suggested the existence of previously unreported 
resonances at about 9 MeV and 23 MeV excitation energy, but their anal- 
ysis did not conclusively identify the multipolarities. 

The work described in this thesis completes the study of 1970y 
discussed above and extends the excitation energy range to 40 MeV. 

Seven giant multipole resonances have been observed. Besides the pre- 
viously reported 7.3, 10.8, and 14.0 states [Refs. 1-4, 18 and 19], the 
existence of resonances at 9.2 and 22.5 MeV are required to explain the 
observed electron spectra. 

Experimental values of the inelastic scattering form factors were 
measured at scattering angles of 60°, 75°, 90°, 105°, and 120° with 
incident electrons of 90 MeV energy. Transition multipolarities have 
been deduced and reduced transition probabilities have been extracted. 
Newly identified are resonances at excitation energies of 9.2 MeV 
(isoscalar EO), 22.5 MeV (isovector E2), and 33.5 MeV (isovector EO 
or E2). Previously reported multipolarities for the 10.8 MeV (isoscalar 
E2) and 14.0 MeV (isovector El) [Ref. 18] resonances have been confirmed. 
The analysis contained here is unable to uniquely assign transition 


multipolarities to the 7.3 MeV (magnetic) and 18.0 MeV structures. 
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II. THEORY 


A. ELECTRON SCATTERING 
The Mott cross section for an extremely relativistic electron 


which scatters elastically from a nucleus of point charge Ze is given 





by [Ref. 6] 
2 z Ss 
aes i ( —- Cos 2 if 
du E; Sin 2 te 


Mott 


The electron is assumed to have spin while the nucleus has neither 
spin nor magmetic moment. Nuclear charge structure in a nucleus of 


finite extent modifies the scattering cross section to [Ref. 6] 


ii-2 





ee 
(se) = (35) |{ ee! ay 
d2 / - 
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Measovored 


e @) is the charge density within the nucleus as a function of T, the 
radius vector from the center of the nucleus, and @ is the momentun 
transfer vector. The Born approximation has been assumed. The term 


describing the nuclear structure is called the square of the form factor, 


( = A 
dw Moeoosvured 


| Fog | : Ge d Ties 


For inelastic scattering, the cross section may also be written as a 
modification to the Mott scattering cross section, but the form is 
more complex, with both transverse and longitudinal contributions 
which are described elsewhere {[Ref. 7]. For inelastic longitudinal 


scattering the form factor is related to the matrix element of the charge 
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operator evaluated between initial and final nuclear states. Although 
transverse components require modification of the simple expression 
II-2, it is conventional to call the experimental form factor squared 
the ratio of measured to Mott cross sections as in II-3. This value 
is compared to theoretical calculations based on particular nuclear 
models. 

For inelastic scattering, the Born approximation may be used to 
compute the theoretical form factor based on an assumed nuclear transi- 
tion charge distribution. First attempts at a solution were made by 
describing the perturbing fields in terms of a multipole expansion, 
and assuming the incident and scattered electron wave functions to be 
plane waves. This so called plane wave Born approximation (PWBA) is 
satisfactory for light nuclei, but becomes less accurate as the number 
of protons in the nucleus increases. For heavy nuclei, the Coulomb 
field surrounding the nucleus distorts both the incident and scattered 
electron wave functions sufficiently that a plane wave solution no 
longer gives the correct results. Consequently, the plane wave solutions 
of the PWBA have to be replaced by phase shifted spherical waves. This 
method, termed the distorted wave Born approximation (DWBA), in practice 
requires computer computations to extract the form factors ([Ref. 8]. 
Form factors thus calculated are plotted as a function of scattering 
angle for each desired electric and magnetic transition. Comparison 
of these plots with experimental results allows for the assignment of 


transition modes to the observed resonances. 


B. GIANT RESONANCES 
The term "giant resonance" has been used until recently to describe 


the giant electric dipole resonance (GDR). This resonance appears as a 
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well known feature in the nuclear excitation spectra produced by such 
photonuclear processes as photodisintegration, photon absorption, and 
inverse reactions such as the (p,y¥ ) reaction. The presence of the GDR 
is characterized by a giant bump in photonuclear spectra at excitation 
energies from 10 to 25 MeV. It occurs throughout the periodic table 
and, according to the Goldhaber-Teller and to the Steinwedel-Jensen 
model, appears at excitation energies that are proportional to the 
inverse nuclear radius, EY = 80A ads (Ref. 9]. 

This conspicuous resonance appears also in inelastic electron 
scattering spectra. The advantage of electron scattering experiments 
over the photonuclear processes for exploring the nature of the giant 
resonances lies in the variable momentum transfer, q. In the case of 
the photonuclear processes, aie momentum transfer is determined uniquely 
by the excitation energy. The giant resonance in the electron scattering 
cross section can be measured at different values of q depending on the 


scattering angle,@, and the incident electron energy, E;; the dependence 


i>? 
being given by 
al all — 
7 = 4&5, 16 alae ae aa 
where E¢ is the final electron energy. In addition, at values of q 
higher than is available from photonuclear work, excitation modes 
other than El become observable. Consequently, the term "giant 
resonance" now refers to all multipole excitation modes which exhibit 
strongly coherent properties. 
Early attempts to explain the giant resonances were based mainly 
on collective models. The initial Goldhaber-Teller model assumed a 
rigid displacement of the protons as a whole against the neutrons 


as a whole, thus producing a large collective dipole moment. The 
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Generalized Goldhaber-Teller model assumes that nuclear matter consists 
of four interpenetrating fluids: neutrons with spin up (nf), protons 
with spin up (pt), neutrons with spin down (n}), and protons with mes 
down (p¥). Any two of these fluids may oscillate 180° out of phase 
with the other two, giving rise to four possible collective modes. For 
example, nucleons with spin up could oscillate against nucleons with 
spin down producing a magnetic resonance ([Ref. 5]. 

The Steinwedel-Jensen model or hydrodynamic model is based on the 
oscillation of neutron and proton fluids within a rigid nuclear surface. 
An extension of this model is the dynamic collective model [Ref. 9]. 
The dynamic collective model replaces the rigid nuclear surface of the 
Steinwedel-Jensen model by a surface that may vibrate. These surface 
vibrations, however, are much slower than the internal giant dipole 
oscillations so that the approximation can be made that the dipole 
oscillations take place in a deformed nucleus. 

A different model of the dipole giant resonances,used mainly for 
light, closed-shell nuclei, is based on the shell model of the nucleus. 
This particle-hole (ph) or Brown model arises from the observation by 
Wilkinson [Ref. 10] that a lp nucleon, promoted into the 2s-ld shell, 
will produce dipole strengths large enough to explain the observed giant 
resonance cross sections. The shell model description [Ref. 10] gave 
the strength of the absorption correctly, but the resonance energy was 
incorrect, and the effect of the particle-hole interaction was noted 
by Brown [Ref. 11]. 

It has been shown [Ref. 9] that such a shell model is separable 
in the relative co-ordinate of all the protons versus all the neutrons; 
and that, upon absorption of a dipole photon, a particular linear com- 


bination of single-particle excitations is produced which corresponds 
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to a motion of all the protons against all the neutrons; i.e., a 
Goldhaber-Teller collective vibration. Hence, the shel] model and col- 
lective models are equivalent in their description of collective states. 

Giants resonances, then, may be described as strongly collective 
modes of excitation in which a considerable number of nucleons take 
part. In the last several years, much progress has been made in the 
investigation of these states using models such as those described 
above. Several modes of excitation other than the GDR have been identi- 
fied, including the giant electric quadrupole resonance (GQR) and the 
giant magnetic dipole resonance (GMDR). A short summary of some of 
these works follow. 

Early electron scattering investigation of the giant resonance 
structure was carried out by Isabelle and Bishop on 16) TRetea ale 
By studying the angular dependence of cross sections, the giant dipole 
was observed. These experiments demonstrated the feasibility of systematic 
Studies of the level schemes and transitions in light nuclei. 

Yamaguchi et. al. [Ref. 12] have carried out extensive investigations 


if oO By making 


of the fine structure in the giant resonance region o 
measurements at both backward and forward angles while keeping the 
momentum transfer constant, they were able to extract the longitudinal 
and transverse form factors as functions of excitation energy. The 
resolution of the fine structure in this manner revealed seven peaks 
associated with the longitudinal form factors and nine with the trans- 
verse. The dipole strength is distributed among a number of states to 
produce this fine structure in the giant resonance region. 


An investigation of the isotopic characteristics of the giant dip ie 


resonance has been carried out by Carlos et. al. [Ref. 13], The giant 
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dipole resonances for the Neodymium isotopes LA? a to gis cs - 


SN, and 
150Nq were studied with a monochromatic photon beam by measuring the 
partial photonuclear cross sections o(¥,n) o(¥, pn) andw(¥,2n). 
As the number of neutrons in the nucleus increased, the dipole resonance 
began to broaden. For sy a the resonance actually split into two 
distinct peaks. An interpretation is given by the dynamic collective 
model in which there is a coupling between dipole oscillations and 


1504 characterizes 


surface vibrations. The splitting of the resonance in 
a deformed nucleus, and the experimental cross section was closely 
reproduced by superposition of two separate Lorentz lines whose centers 
are separated by four MeV. In the framework of the hydrodynamic model, 
one can assume for an axially deformed nucleus two modes of excitation 
at differing energies, equivalent to oscillations in the direction of 
the two axes. 

Pitthan and Walcher [Ref. 14] did (e,e') investigations with Ce, La, 
and Pr. All three elements show distinct resonances at excitation 
energies of 9, 12, and 15 MeV. The resonance at 15 MeV was identified 
as the El giant resonance based on correlation of the excitation energy 
and strength with the GDR seen in photo-neutron experiments. The 12 
MeV resonance exhibited longitudinal properties similar to the El 
but did not appear in (@,n) studies. It was assigned an E2 or EO 
excitation mode based on comparison to form factors developed from the 
hydrodynamic model. The 9 MeV resonance was found to be transverse 
and most likely an Ml transition. 

Later, Lewis and Bertrand [Ref. 15] observed that the (p,p') spectra 


27 54 1205, | aa 209 


et -Al, Fe, Bi consistently produced enhancements at 


excitation energies slightly lower than those of the known giant 
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dipole resonance. Since the resonance energies of the giant dipole had 
been previously well measured, they concluded that it was unlikely that 
the giant dipole excitation was responsible for most of the observed 
enhancement. This apparent resonance shift was due to the presence of 
a giant isoscaler-quadrupole vibration located at EY. = 63 a /3 MeN. 
In forward angle (e,e') experiments, Fukuda and Torizuka observed 
{Ref. 16] similar El and E2 ragenaaces in 900. The usual dipole reso- 
nance was seen at 16.65 MeV and additional resonances were seen at 14 MeV 
and 28 MeV. Based on a prediction of the hydrodynamic model that a 
quadrupole giant resonance should appear at a position 1.6 times the 
peak energy of the usual giant dipole resonance [Ref. 17], the 28 MeV 
resonance was assigned an E2 transition. An assignment of either EO 
or E2 was made for the 14 MeV peak. Based on similar results for Soke 


116 208 


Sn, and Pb [Ref. 16] these new giant resonances seem to be universal 


features existing at excitation energies of 65 AV 3yey and 130 AV 3yey 


1/3 dependence of 


in all nuclei. Table I contains a listing of the A” 
the resonances identified in this work and a list of references where 
others have seen or predicted corresponding resonances in 197 ay and 
other nuclei. 

TABLE I 
=/5 


Dependence of Yk Resonances 


A 





References 







t-[Refs2 2b°G 22) - Ce,aLa, Pr [Ret 25) 





Au [Refs. 3 & 18] 
Au [Reis 3,-18,. and: 25] 
Pb [Refs. 24 § 26] 
Au [Ref. 3] Pb [Ref. 24] 






















t [Ref..22) Fb [Ref. 26] 


t = Theoretical Prediction 
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C. TRANSITION PROBABILITIES AND SUM RULES 


The transverse-electric and eee form factors are given by [Ref. 27] 


2 1) poe Blt. bay, 7d, ) II-5 





2 
| Fe ue Cad] = = 


The continuity equation provides similar relationships for the longitudinal 
form factor. The coefficients B are the reduced nuclear transition prob- 
abilities and, as such, contain all of the nuclear physics of the inelastic 
scattering problem. The theoretical form factors are calculated by the 
GBROW computer program as described by Ziegler ([Ref. 27]. 

To determine whether or not the resonances observed are collective 
phenomena, comparison of experimental to single-particle reduced tran- 
sition probabilities may be made. For a heavy nucleus, one would expect 
a strength greater than several times the single-particle value. Further- 
more, a giant resonance should exhaust an appreciable fraction of the 
appropriate sum rule. If an observed resonance greatly exceeds the sum 
rule under the assumption of a particular transition multipolarity, it 
is unlikely that the multipolarity assignment is correct. 

The single-particle reduced transition strengths are given by the 


Weisskopf unit [Ref. 28], 
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-1/3 197 : 
where R= 1.2 A = 6.98 fm for Au. Table II provides values of 


the single-particle transition strengths for Tt hae 

Another, to some extent model-independent, evaluation of the transition 
strength of EL modes may be achieved by expressing the strength relative 
to the energy-weighted sum rule (EWSR). Sum rules for magnetic transitions 
are more dependent on the structure of the nucleus in question and, there- 
fore, no comparisons to these have been made here. States at excitation 
energies less than 15 MeV appear to be AT = 0 except for the EL mode 
which is AT = 1. The sum rule for isoscalar (AT = 0) excitation modes 
with L?1 is given by [Ref. 29] 


z 2b-2 
ee ve Seley ace > 
S (EL 4T=0) = 215 BCR Es Saar . I1-8 


The isovector (AT = 1) sum rule for L>1 is related to the isoscalar 


sum rule by 


N 
S (EL T=!) S (EL \OT=0) es Tks 


The corresponding sum rule for an isoscalar monopole (E0) excitation is 


[Ref. 30] 


S (EO) = ete | Me;{ 


lt 


+2 z 
SLA is ON | II-10 


where Mes is the monopole matrix element. The strength of El resonances 
can be expressed relative to the sum rule given by Warburton [Ref. 31]: 


ch ( ay bbe 
Seer 2 ie tan A; II-11 





197 : oe 
Table III gives values of these sums for Au using weighted radii of 


the ground state charge distribution. 
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TABLE II 


Weisskopf Units for “ig 


6.6 


3.35x102 


1.63x10° 1.44x102 


* Units are fm2h 





TABLE III 


Sum of EWSR For 197 ay 


9.2x104 me een 


7.04x102 e--MeV-fimz 
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6.15x10° e-Mev-£m® 


9.16x10° e--Mev-fm® 
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D. QUASI-ELASTIC SCATTERING 

Electrons incident upon a nucleus may scatter elastically with 
individual nucleons within the nucleus. This process is called quasi- 
elastic scattering. If sufficient momentum is imparted to the nucleon, 
it will be knocked out of the nucleus with some kinetic energy. The 
simultaneous energy loss by the incident electron may result in a 
measurable quasi-elastic peak. It is of interest to see at which value 
of EY this contribution appears on the spectrum of a= , to see if 
it contributes to, and therefore must be subtracted from, the scattered 
electron spectrum in the region of the giant resonance. 

The inelastic electron-nucleus cross section has been calculated 
in three energy-loss regions with the Fermi gas model for the nucleus 
[Ref. 32]: the quasi-elastic region; the threshold pion region with its 
pion electroproduction cross sections; and the region containing the 
3-3 resonance of the proton. Pion production begins at 135 MeV and the 
3-3 resonance is 298 MeV above the elastic peak, so these regions are 
not of interest in this work. 

In the Fermi gas model, the electron scatters elastically from a 
single nucleon in the free Fermi sea, with the recoiling nucleon required 
to scatter out of the Fermi sphere. Because of the Pauli principle, the 
nucleon cannot scatter into an already occupied state. The final nucleon 


energy is given by 
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where k is the nucleon momentum and M is the nucleon mass. The initial 


(bound) nucleon energy is 


os ie Ui eS Be + Ule) 5 P1213 
Kt 2M 2M 


where the effective nucleon mass, M, is given by [Ref. 33] 


* 


MN ze os II-14 
ies 

and U(k2) is the effective single-particle potential in nuclear matter. 
This potential effectively shifts the electron energy loss, WJ , to take 
into account the nuclear binding. Moni treats U(k?) as a constant, E * 
equivalent to the average nucleon interaction energy. The quasi-elastic 
peak energy loss will increase with q and approach the effective free 
nucleon kinetic energy, an / (2M") , at large values of q. 

The magnitude of the above shift for the gold nucleus is of interest. 
The average nucleon interaction energy, E , and the nuclear Fermi 
momentum, k. , have been determined for nine elements from C14 to 
208pp [Ref. 34] by a least squares fit of theory to quasi-elastic peak 
data. In the cross-section formula, an energy-conserving delta function 


appears, given by 


Fi a ea) oe SD 
SO om = ae se). II-15 
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which determines L) . The nucleon momentum is approximated by the Fermi 
momentum and is roughly constant at 260-265 MeV/c for heavy nuclei from 
nickel through lead. E i a 43 MeV [Ref. 34]. The energy-conserving 


delta function requires that 


hed ae 
tehoen tt: innards yaveFn 2 II-16 
NM 2M 
For the smallest momentum transfer of these experiments, W = 73 MeV. 


The threshold value is G)= jose 43 MeV. Thus explicit corrections 
for the quasi-elastic contribution need not be made for the spectra of 


this work. 
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III. DATA ANALYSIS 


A. DATA REDUCTION 
Five experiments were performed with an incident beam energy of 


©; and 120°. 


90 MeV and at scattering angles of 60°, 75°, 90°, 105 
The experimental information is given in Table IV. In general the data 
analysis procedure described by Warshawsky and Webber [Ref. 1] was 
followed throughout this work. The “ee choice of resonance energies 
was determined from known resonances [Refs. 1, 2, 3, and 4] and by visuai 
inspection. Previously determined resonances occur at 756° Oc25 12.0, 
and 14.0 MeV. Visual inspection of the background subtracted data 
revealed wide structures at 18.0, 22.5, and 33.5 MeV. 

The fitted inelastic spectra extended from 4 to 40 MeV, except 
the 120° spectrum which extends only to 31.5 MeV. In some cases the 
data extends to 50 MeV but no structures were observed at this high 
excitation energy. Since the 120° data has different boundary conditions 
for matching the background at high excitation energies, only cross 
sections for resonances below 20 MeV were extracted at this angle. 

Breit-Wigner resonance forms were assumed for each resonance and 
computer code NAW (see appendix) simultaneously fits the entire inelastic 
spectrum including the radiation tail and background. The new resonance 
positions and widths were manually varied while the NAW code varied the 
peak strengths and background parameters until a best fit was obtained. 

The criteria used to determine a reasonable fit were as follows: 

1. The data and calculated spectrum should coincide visually. 

Zs < per degree of freedom should be less than one. The data is 


not strictly statistical because the detector momentum interval is larger 
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than the momentum increment of the epoctronceer field and hence correla- 
tions exist between energy bins. 

3. The difference between the total b Wes for successive computer 
iterations was aX? ¥ 0.5. Typically there were 300 degrees of freedom 
and a change of 0.5 in the total chi squared resulted in 17 x 19-2 per 
degree of freedom. 

4. All observed resonances should consistently fit each of the 
five spectra. Table V lists the positions and widths of the resonances 
required. 

Figures 1 - 5 present the experimental data corrected for spectrometer 
dispersion effects with the fitted total background and the individual 
resonances drawn. Figures 6 - 10 present these spectra with the total 
background subtracted, revealing the resonance structures of Table V. 

The areaSunder the inelastic resonances were calculated from the 
Breit-Wigner resonance shapes and the ratio of each inelastic resonance 
area to the area under the elastic peak was calculated (A;/A,). The 
inelastic form factors squared are obtained by multiplying A./Ke by 
the square of the elastic scattering form factors Gar Table VII 
presents the experimental ratios of the inelastic to elastic areas and 
the square of the inelastic form factors. The elastic form factors 
were calculated using the program of Rawitscher and Fischer [Ref. 35]. 


The charge distribution parameters were c = 6.38 fm and t = 2.32 fm. 
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TABLE IV 


Experimental Conditions 


Target 
Thickness 


.048 g/cm? 


.096 g/cm” 


.096 g/cm? 
£144 g/cm? 


.240 g/om? 


TABLE V 


Incident 
Energy 


90.32 MeV 
90.24 MeV 
90.21 MeV 
89.98 MeV 


89.60 MeV 





Observed Resonances 
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TABLE VI 


Scale Factors for Converting 
Arbitrary Units to Experimental Units 


Scale Factor MeV-gr / Arb. units 


9.9 x1074 


1.89x1074 


2.99x107> 


1.53x107° 


1.47x10° > 
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TABLE VII 


Area Ratios and Inelastic Form Factors Squared 
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TABLE VII Continued 
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B. ERROR ANALYSIS 

The procedure outlined in the previous section provides values of 
the ratios of Aj/A, but the uncertainties attached to these ratios 
requires additional considerations. The program NAW computes only 
statistical errors. There are other errors, both systematic and random 
which are not strictly statistical. Some of these are the errors result- 
ing from the instrumental effects and are outlined by Bernard and 
Praverso [Ret. 36]. 

To obtain a more realistic estimate of the errors in the ratios 
A,/A,: the non-statistical uncertainties are assumed to be those which 
arise from the uncertainties in locating the resonance positions and 
determining their respective widths. These arise from two sources: 
the overall instrumental uncertainty which is estimated to be 0.15 MeV, 
and the freedom allowed in the resonance peak positions and widths in 
fitting these spectra. The total uncertainty of the energy in this work 
is the sum of the fitting possibilities and instrumental errors. 

The Breit-Wigner shape is then analysed to see how these variations 
in energy effect the area under the resonances. The Breit-Wigner shape 


is given by 


Ao (EEM)  _ A 
dE (EE ier ee 
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where ER is the resonance energy, ‘e is the resonance width, A is the 
strength parameter, and E is the excitation energy. The resonance cross 
section is given by the integral of the Breit-Wigner form 


-|, Ep 
eter) | ae i tus kcal) 





EII-zZ 


The variations in position and width introduce an uncertainty in 


the resonance cross section, 


oo 
Ar(e.?)= [Se 





La 
E A, : 
AE, ee III-3 


The absolute values of the partial derivatives are assumed since the 
peak position and width uncertainties may be correlated. 

The fractional error in the resonance a section is 

2 Tv 2 Tan Ct) oo gp BE pe 2a \ 

2 [ree |4* “ [ se ae Ee oes : 

SS E47 
A i see 

¥ iV 
This uncertainty is then added to the statistical error to provide 


the total error. Table VIII lists the uncertainties from the two sources. 
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TABLE VIII 


Fractional Error Contributions 
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Stastical Errors 
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IV. DISCUSSION 


me | COLLECTIVITY 

Data analysis techniques explained in Chapter III were used to sub- 
tract the background and extract the elastic and inelastic spectra 
. of the scattered electrons. A least squares fit was used to evaluate 
these spectra and resulted in the identification of seven giant 
Fesauances at excitation energies of 7.3, 9.2, 10.8, 14.0, 18.0, 
22.5, and 33.5 MeV. These resonances have strengths corresponding to an 
appreciable fraction of the appropriate sum rule as exhibited in Table IX. 
Additional support to their collective nature is given in Table X where 
it is seen that the transition strengths are greater than several Weis- 


skopf units. True giant resonances are therefore being observed. 


B. MULTIPOLARITY ASSIGNMENTS 
7.3 MeV Resonance 

The multipolarity assignment of the resonance at 7.3 MeV is un- 
certain. Buskirk et. al. [Ref. 18] and Lone et. al. [Ref. 19] have 
observed this structure in backward angle (e,e') experiments and indicate 
a transverse excitation. A comparison of the data with calculations using 
GBROW [Ref. 27] and a code by Drechsel [Ref. 37] is incapable of deter- 
mining a definate multipolarity assignment. Among possible reasons for 
this situation are a multipolarity admixture, existence of unknown 
longitudinal resonances influencing the forward angle data, and experi- 
mental difficulties encountered by the high counting rate at forward 


scattering angles. 
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9.2 MeV Resonance 

The analysis of the 9.2 MeV resonance (53 a3) favors an EO 
multipolarity assignment. This resonance exhausts 38% and 18% respect- 
ively of the EO and E2 energy-weighted sum rules and its strength 
corresponds to five Weisskopf units when an E2 assignment is made. 
However, the angular dependence of the experimental form factors favors 
EO (Fig. 12). The result would be in agreement with the isoscalar 
giant monopole resonance predicted at 56 AV Btey {[Ref. 21] and at 
58 v/s MeV [Ref. 22]. Resonances exhibiting the characteristics of 
EQ or E2 transitions were also seen at Ex = 53 AV Ayoy an N= 82 


ae [Ref. 24]. An El multipolarity assignment 


nuclei [Ref. 23] and in 
is ruled out by comparison with DWBA form factor calculations and further- 
more by the fact that no resonance has been seen in 0B - absorption 
measurements [Ref. 2] at this excitation energy. 
10.8 MeV Resonance 

The resonance at 10.8 MeV (63 a i/3, was seen previously in 197 Ay 
[Refs. 3 and 18] and in many other nuclei (see e.g. [Ref. 5])}. Generally 
an assignment of E2 was favored in these cases. In this experiment a 
transition strength of 14 Weisskopf units is obtained assuming an 
electric quadrupole excitation. Comparison of the experimental form 
factors to the angular distribution calculated from the DWBA form factors 
is inconclusive in distinguishing between an EO and E2 modes (Fig. 13). 
However, the experimental strength corresponds to 113% and 68% of the 
EO and E2 sum rules, respectively. Together with our results on the 
9.2 MeV resonance, the E2 assignment is therefore favored. This resonance 


is thought to be the isoscalar quadrupole state predicted by the 


collective vibration model of Bohr and Mottelson [Ref. 38]. 
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14.0 MeV Resonance 
The 14 MeV (81 a7l/3) giant dipole resonance was measured with 
monochromatic 2&- experiments [Ref. 25] and has also been seen in other 
(e,e') experiments [Refs. 3 and 18]. For this experiment, form factor 
comparisons to DWBA calculations show reasonable agreement at forward 
angles, but exhibit deviations at backward angles. An attempt was made 
to further explore the nature of the giant dipole transition by assuming 
both volume and surface oscillations separately in the DWBA calculations 
(Fig. 14) [Ref. 23]. Transition strengths of six and eight Weisskopf 
units and 84% and 104% of the sum rule are obtained for the volume and 
surface oscillation modes respectively. However, neither model adequately 
describes the observed data. As suggested by Pitthan [Rei 29) sine 
discrepancy at backward anglesmight be explained by an electric spin flip 
contribution to the El resonance [Ref. 40]. Using only the forward angle 
data to extract the reduced matrix elements, the volume oscillation 
strength exhausts the sum rule (98%), whereas the surface oscillation 
exceeds the sum rule. Thus a volume oscillation is favored. 
18.0 MeV Resonance 
A resonance of undetermined multipolarity is observed at 18.0 MeV 


-1/3 


(105 A ) excitation energy. Nagao and Torizuka [Ref. 24] have observed 


een but the state in 197 Ay 


E3 strength at the corresponding energy in 
does not conform to this multipolarity. (0-3) and M (1-2) calculations 
‘have been compared to the data (Fig. 15), but no reasonable agreement 


is found. It appears that the cross-section at this excitation energy 


contains a mixture of various transitions. 
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22.5 MeV Resonance 
Comparison of the experimental results for the 22.5 MeV (130A7 1/3) 
resonance with DWBA calculations (Fig. 16) gives support to an E2 
assipnment. At forward angles the data is indistinquishable between 
the EO and E2 transitions, but the 105° measurement is incompatable 
with the EO assignment. The Sendai group has reported [Ref. 3] an 


197 A and 


electric monopole or quadrupole at this excitation energy in 
also at the corresponding energy in 208pp [Ref. 24]. Other multipolarity 
assignments are unsuccessful in predicting an angular distribution in 
agreement with the data. The 22.5 MeV state is thought to be the 
isovector quadrupole state predicted by the collective vibration model 
of Bohr and Mottelson [Ref. 38]. 
33.5 MeV Resonance 

A new resonance, previously unreported in tt hka, is observed at 
33.5 MeV (194 Avl/3y, This corresponds roughly to the isovector giant 
monopole resonance predicted by Suzuki [Ref. 22] through calculations 
based on the collective vibration model of Bohr and Mottelson. Both 
the EO and E2 assignments exceed their respective energy-weighted 
sum rule prediction. However, our line shape fit, while very sensitive 
to any structure in the spectrum, may tend to over-estimate cross- 
sections for very wide resonances. For this reason, the B-values may be 
too large. Comparison of the experimental form factors to the DWBA 
“calculations allow both EO and E2 multipolarities to be possible. A 


208 : 
Pb at a corresponding energy of 


Similar new resonance was observed in 
32.5 MeV [Ref. 26] and tentatively given an EO or E2 assignment. It is 
to be emphasized that consistent fits for all the spectra 97 au) 

over the entire range from 5 MeV to 40 MeV excitation energy were only 
possible with the inclusion of this resonance at 33.5 MeV. 
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Transitions Strengths in Weisskopf Units 
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V. CONCLUSIONS 


Giant multipole resonances in 197 ny were studied with inelastic 
electron scattering. Resonances at 7.3, 9.2, 10.8, 14.0, 18.0, 22.5, 
and 33.5 MeV excitation energy were observed. Of these, the 9.2, 18.0, 
22.5, and 33.5 MeV states have not previously been reported in the open 
literature. Multipolarity.assignments of these states were made as 
was discussed in section IV-B. 

The 14.0 MeV state is identified as the much studied giant dipole 
resonance, see e.g. [Refs. 3, 5, § 25]. Our analysis suggests that this 
El state may be described by a volume oscillation of the charge distribution. 
Furthermore, the data is compatible with a suspected transverse spin flip 
contribution to the El cross section. 

The angular distributions of the cross-sections for the 10.8 MeV 
and 22.5 MeV resonances are compatible with previous E2 multipolarity 
assignments [Ref. 3]. In addition, this investigation adds evidence 
for the exclusion of a monopole assignment for the 10.8 MeV state. 

A state was found at 18 MeV in mc In 208p, a resonance was 
seen [Ref. 26] and cross-section analysis of another (e,e') experiment 
found E3 strength between 16 MeV and 22 MeV [Ref. 24]. The analysis 
of the experimental results for this energy region in 197 ay is not 
capable of assigning a definite transition multipolarity. The data 
suggests admixtures of longitudinal and transverse contributions to the 
cross-section. 

The 7.3 MeV transition is known to be magnetic. However, this 
investigation, like others [Ref. 19], was not able to determine the 


transition mode. 
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Theoretical Paes een concerning tie existence and location of 
the giant monopole resonance is found in the literature [Refs. 21 and 22]. 
New resonances, discovered in this experiment at 9.2 MeV and 33.5 MeV, 
have angular distributions in agreement with monopole assignments. Their 
excitation energies are comparable to theoretical predictions [Ref. 22] 
of the location of the isoscalar and isovector monopole states. For 
the lower state, an EO transition is favored over an E2 assignment. 
However, for the 33.5 MeV state, the transition multipolarity is 
inconclusive. The excitation energies and transition multipolarities 


are presented in Table XL. 
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TABLE XI 


Multipolarity Assignments 


Excitation Transition 
Energy (MeV) Multipolarity Comments 


undetermined probable magnetic, 
mixture 


EO Isoscalar 
E2 Isoscalar 


El Isovector, volume 
oscillation 


undetermined 


E2 Isovector 


EO or E2 Isovector 
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APPENDIX A 
PROGRAM NAW, SPECTRUM FITTING CODE 


A. GENERAL FEATURES 

Naw is a Fortran program developed at Darmstadt [Refs. 41 § 43] 
which uses analytical expressions to approximate the inelastic and 
elastic line shapes and radiation tail to best fit the experimental 
spectra. The elastic line shape is assumed to be two half gaussians of 
differing half widths matched to a hyperbolic shape approximating the 
radiation tail. Work done by the experimenters at the Institut Fir 
Technische Kernphysik at Darmstadt have found this parameterization to 
be an adequate estimate of the shape of the elastic peak (see Figure 18) 


[Ref , 41]. 
gaussian 


gaussian eh 


hyperbolic 





High E Low E 


relative zero 
peak position 


ASSUMED ELASTIC SPECTRUM SHAPE 
FIGURE 18 


The elastic radiation tail results from processes occurring during 
elastic scattering that produce lower energy electrons. The radiation 


tail is calculated from the theory [Ref. 43] at several designated energies 
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and these theoritical values are fitted with a polynomial of the fourth 
degree. The total background, which includes the radiation tail, is then 
simultaneously fit with a seivits of Hreiewiner line shapes representing 
the nuclear excitation structure to reproduce the experimental spectrum. 
The positions and widths of the resonances are free parameters to be 
adjusted by best fit or to be assigned from known values. The areas 
under these resonances are then calculated from which the ratios A;/Ag 
(inelastic peak area divided by the elastic peak area) are determined. 
The values A;/Ag are then multiplied by the calculated scattering cross 
section resulting in the experimental inelastic cross section for each 


of the resonances. 


B. CONTROL CARDS 

NAW uses the output from program KIK as the input data. Since 
NAW's first use at che NPS IBM 360/67 [Ref.1] it has undergone various 
modifications. The present form is reproduced on the following pages, 
along with the control information. NAW accepts 200 elastic peak data 


points, 500 inelastic data points and 25 inelastic resonances per 


spectrum. 
NAW CONTROL 
A TEXT Alphanumeric identifying remarks (20A4) 
2. NORM, PN,, PNz, PNy, PNg, PN» PN (7F10.5) 
NORM Spectrometer energy normalization factor to change the 


origional energy units (OEU) to MeV: [MeV] = NORM*[OEU] 


PN; ie2-7, Starting values for the various parameters being 
used to fit the elastic spectrum with an analytical 
expression. 

PN, Half-width of a paussian to fit the high energy side 


of the elastic peak. 
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3 Half-width of a gaussian to fit the low energy side of 
the elastic peak. 


PN, Matching point for the low energy gaussian and the 


hyperbolic segment. PN, is specified in MeV from the 
position of the zero point and is a negative number. 


PNs Dummy . 


PN¢ PN 0, Height of a second elastic resonance peak if 
necessary. 

PNo Position of a second elastic resonance peak if necessary, 
only with PNg* 0. 

EE, THETA, D, UNT, (4F10.5) ~ 

EE Energy of the elastic peak in MeV. EE is used as a 


starting value in the elastic fit. 


THETA Electron scattering angle in degrees. 


D Target thickness in g/cm. 

UNT Constant background at the high energy end of the elastic 
resonance spectrum in counts/ac. 

EEL,,, EEL, GSSW , GSSW,, DSSW, Eo¢ (6F10.5) 

EEL, High energy limit for fitting the elastic peak in OEU. 

EEL, Low energy limit for fitting the elastic peak in OEU. 


GSSW | Upper limit for calculating the radiation tail in OEU. 


GSSW, Lower limit for calculating the radiation tail in OEU. 


DSSW Energy step to be used in calculating the radiation 
tail in OEU. 


E ef See card 7, Lp. If Evef = 0, Eref is calculated internally. 
aw Fl (2F10.5) 
Ey High electron energy limit for fitting the inelastic 


spectrum in excitation energy (MeV). 


Ej Low electron energy limit for fitting the inelastic 
spectrum in excitation energy in (MeV). 
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6. GR, GR, or2, cR2, cRé, cr? (6F10.5) 
Regions where it is expected that only background 
contributes to the spectrum; i.e., between resonances 
in the inelastic spectrum. There must be at least one 
set of GR?,, GR} input, with a maximum of three. It is 
recomended to set GRyxGSSWe GRimxGSSWu for giant 
resonances, energy in MeV. 


1.) Le -2=l to 22 (2213) 


L: 1 A continuous curve consisting of two gaussians and a 
segment of a hyperbolic is used for both elastic and 
inelastic fits. 

2 The elastic resonance will be fit as if L,=1 and the 
inelastic resonance will be fit using the Breit- 
Wigner formula. Each curve will be multiplied by a 
straight line with an equation Y=1 +P, (E,-Ex)). oe 
will be Sees on card #10, Ex is the excitation 
energy, and E, is the position of the maximun 
as specified on card #8. 

3 The elastic resonance will be fit as if L}=2 except 
that NAW will determine the values P. and card #10 
will be left blank. : 


Lo Parameters which choose the analytical expression to 
fit the background due to the radiative and experimental 
processes, L2 is a two digit number, the units give the 
order of a polynomial e.g. unit=3, P,+PjA+P,A°, the P;'s 
are free fitted paramters. The decades chodse the 
function to be multiplied and the argument to be used. 


2 - ial 
eee eee een 
10< L,< 20 P,+PoA+P A+... FSSW 


2 ef ] ae 
20 < Ly < 30 {Py+P2A+P 3A a eee Fsswe exp Ce | 9 | 
2 y I.E 
30< L,< 40 {Py +P AHP SA +...4 FSSW f exp [Pizs ey 


FSSW is the calculated radiative tail, E is the energy 
of the outgoing electron, E., is the energy of the 
elastic peak. A-(E-E,.-) for L,= 50, EL oe is read in 
oncard $7. If SO<L =< 90 the“analytical form correspond- 


ing to L5=L,-50 is chéosen, but A is replaced by 1/A. 


1¢ L,< 10” Py TPAtP 


L The number of resonances with fixed energies. If La< 0 
then the resonance will be kept fixed relative to ELC), 
but EQ) will be fitted. See card #8. 1 


Ly The number of resonances with energies to be fitted. 
Note that |L.| + L, must not exceed 25. 
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The number of lines widths to be fitted. 70 E.g. =3, 
means that three line widths are free and in fact are 
E (1), E.(2), and E_(3). If there are more than three 
résonancés , E,.(3) tiirough E,(N) will all be fitted 
with the same line widths. 


For L, = 2 © .thé first Le peaks will be fitted with 
Breit-Wigner forms, the remainders with gaussian 
shapes. For L.=0, L¢= [La[+ Ly 


Target geometry: 
=Q Transmission -normal thickness entered on card #3 


=] Reflection -normal thickness entered on card #3 
=2Z Effective thickness entered on card #3 


Output control: 

=0 Maximum output 

=1 Minimal output-no tables and no intermediate calculations 
=-| Limited output-no tables 


Plotting control: 


=0 Results are not plotted. 
=Z Results are plotted. 


Repetitions of NAW procedure. 


=0 No repetitions are performed. 

>0 In this case it is possible to evaluate the data with 
new parameters. The new parameters will be entered 
on additional cards placed after the control cards. 


Inelastic spectrum fit iteration control. Only 
override to the default options need be entered. 


Number of iterations where all energies are fixed, 
default L,,=2. 


Maximum number of iterations; default is L739. 
Minimum number of iterations; default is L375. 


Number of iterations where width are held constant in 
case % actual >3* theory. 


=0 The radiation tail will be calculated using the 
Ginsber's method [Ref. 43]. 

<0 The radiation tail will be calculated using Schiff's 
method [Ref. 44]. 
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1 “Lig Radiation tail fit control. Only override options need 
6 be entered. 
16 Analytical representation of the theoretically 
calculated radiation tail; default is Lig. 
2 3 4 
c I= 
1 SSW PS xi PA f PA Gs PAA = PLA 
Ps if 2 : 2 
=2 SSW=P, + PA + PA + PAE * PLE 
2 a : 2 3 
=3 SSW=P) + PoA + PLE + PAE + PLE 
= z 2 3 4 
=4 SSW=P 4 + PoE + P2E + PAE + PLE 
a = 2 3 4 
=5 SSW=P; + PoU + P2U" + PyU™ + Pou 
P. will be fitted internally by NAW 
U-E - Eh 
Ee GSSW cP (GSSW,, = GSSW, )/L 7 
<0 the radiation tail parameters P. and By will be 
specified on card #15. . 
L,, The number of parameters for the radiation tail fit 
: (3 < hiya <= 5); default is Lj = 4, 

lig See Lig; default is Lig =-2. 

Lig 20 Error ellipsoid control values will not be calculated. 

easy 

Loy >0 The calculation of the radiation tail will be 
modified by the factors entered on card #16. 

Lo5 =3 Inelastic spectrum results. Card images will be 
written on file 7 and may be stored on the data cell for 
further handling. 

E (1), E.(2), ---, EC) (SF10.5) 
Resonance excitation energies in MeV; enter those 
with variable positions first. i= |L.| + Ly < 2s. 

MElD, V( Ziad. @) (8F10.5) 
Resonance line widths, i= |L,| + L,< 25 


If any V(i) = 0, the width ofthe elastic line will be 


assumed. 
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PO ee tb 2 as «ng PUL) (SF10.5) 


The P. factors when L = 2 on card #7, i = |L5| + L,<25. 


a1s Z, Ag RAD) FIA CT : (4F10.5, 20X,2F10.5) 
Z Atomic number. 
A Mass number. 
XRAD Radiation length in g/cm? . 
FIRA Elastic form factor for E, and THETA of the experiment. 
C,f The fermi charge distribution parameters. Since these 


12. XMAS, 
XMAS 


14. U 


ih. Pb, Pe. 


Pl=Ps 


H 


values are used as the parameters for the trapezoidal 
charge distribution potential used in the radiation tail 
calculations in Born Approximation, these values must 
be entered. 


(2F10.5) 
Nuclear mas equivalent in MeV. 
Isotope percentage (100% = 1.00). 


‘the A. are the coefficients for the calculations of the 
ejastic form factor using a power series expansion on 
qe 
2 ° 2 (i 
= 2 L\(l— 
Fai x Aq? G=1) 
When A,#0 the elastic form factor will be calculated 
using the Born Approximation. 


(8F10.5) 


Ujmay be used to get better starting values for the fit 

by taking the background parameters from a previous fit. 
If the U;=0, the starting values are chosen internally. 

See. card #7. 


Baa eS) be (6F10.5) 


This card is read only if Ligs 0. 


Radiation tail parameters, see L,. for the expression. 
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16. XSSW, -XM, XB, (3F10.5) 


XSSW 


This card is read only if Lj; >0. 

Factor Multiplying the area under the elastic peak, a 
e 

Factor multiplying the Mé@ller tern. 


factor multiplying the B 
The equation for the radiation tail is now given by 


SSW = a _ * Miler term) + XB(Brems. term)* (XSSWA, 
MAXI are fers factors as described by Maximon and Isabelle. 


The data are assigned to file #4. The elastic and the inelastic 


data must be followed by one card image having the data energy step 


width entered beginning in colum 72. Repeated evaluations of the same 


set of data (i.e. L,, = 0) must have the following set of cards placed 
10 


after the 14 to 16 control cards as described above; one set for each 


repetition to be performed. 


1) 
2) 
3) 
4) 
5) 


Ly erate Lig see card #7 
Ey, Eq see card #5 
EO) Weiss EG) see card #8 
=f 6 Pe (gre) see card #9 
P(L) 5 2, bas) see card #10 
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